Disorder and Quasi-Particle Interference in Heavy Fermion Materials 
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Using a large- N approach, we study the effect of disorder in the Kondo-screened phase of heavy- 
fermion materials. We demonstrate that the strong feedback between the hybridization and the 
conduction electron charge density magnifies the effect of disorder, such that already small con- 
centrations of defects strongly disorder the materials' local electronic structure, while only weakly 
affecting their spatially averaged, thermodynamic properties. Finally, we show that the microscopic 
nature of defects can be identified through their characteristic signatures in the hybridization and 
quasi-particle interference spectrum. 



Introduction Heavy fermion materials are characterized 
by the presence of localized degrees of freedom, i. e. mag- 
netic moments residing on the rare earth or actinide ions 
and itinerant electronic states formed by spd-electrons of 
the system and the strong correlations between these de- 
grees of freedom^. The resulting interplay gives raise to a 
wide range of possible ground states ranging from various 
magnetic and superconducting phases^^to semiconduct- 
ing and metallic phases^ with strongly enhanced quasi- 
particle mass, or without well-defined quasiparticleaS and 
even to enigmatic phases, where the order parameter is 
not known despite decades of intense research^. Only 
very recently, scanning tunneling spectroscopy (STS) ex- 
periments have succeeded in probing the local electronic 
structure of several heavy fermion compounds, such as 
URu 2 SiJSHIol, YbRh 2 SiM CeCoIn 5 and CeRhlnrP, at 
sufficiently low energy and temperature to infer ground 
state properties. In particular, by utilizing the spatial os- 
cillations in the differential conductance, dl/dV, induced 
by defects and performing a quasiparticlc-intcrfcrencc 
(QPI) analysis, it was possible to ma p out t he electronic 
band structure near the Fermi energy^- 12 * 20 !. STS experi- 
ments also reported defect-induced spatial oscillations in 
the hybridization possessing a wavelength that is deter- 
mined by the unhybridized, small Fermi surface of the 
conduction band^, thus confirming an earlier prediction 
by two of the authors le . Surprisingly, the same STS ex- 
periments also found that already a small, 1% concen- 
tration of defects strongly disorders the hybridization in 
the entire system, while in general thermodynamic bulk 
measurements are largely insensitive to doping levels up 
to a few percenlP^l Resolving this apparent contradic- 
tion between STS and thermodynamics clearly requires 
a more microscopic understanding of defect-induced ef- 
fects in heavy fermion materials. 

In this article, we address this issue and show that the 
strong feedback effect between the defect-induced spatial 
oscillations in the hybridization and the charge density of 
the conduction band enhances the strength of the disor- 
der, thus explaining the strong disorder effects observed 
experimentally. At the same time we find that thermo- 
dynamic properties of the system, such as the spatially 



averaged specific heat, are only weakly affected by de- 
fect concentrations of a few percent, thus explaining the 
qualitatively different STS and thermodynamic observa- 
tions. Moreover, we show that the spatial form of the 
hybridization oscillations as well as the QPI spectrum 
provide insight into the nature and scattering strength 
of the defects. 

Model To study the effects of defects in heavy fermion 
materials, we consider the Kondo-Heisenberg Hamilto- 
nian 

H=~t ^ C L C r',<x ~ Vj24,a c r,<r 
{r,T>),<T r,cr 

r (r,r') 

where t is the hopping element between nearest-neighbors 
of a two-dimensional square lattice describing the light 
conduction band, and c\ a (c r CT ) creates (destroys) a con- 
duction electron of spin a at site r. We choose a chemi- 
cal potential of fi — — 3.618t, resulting in a Fermi wave- 
length Xp = 10 (we set the lattice spacing, ao, to unity) 
and an electron density of n c ~ 0.062 of the (decou- 
pled) conduction band. J is the Kondo coupling be- 
tween the conduction electrons and the S = 1/2 spin 
of the magnetic atoms , and I is the strength of the 
nearest-neighbors antiferromagnetic interaction between 
the magnetic atoms, which is treated here as an inde- 
pendent coupling constant)^. Moreover, we consider de- 
fects in the form of missing magnetic atoms (i.e., vacan- 
cies), and non-magnetic atoms that are substituted for 
magnetic ones, where the latter lead to potential scat- 
tering term U Yin. a c r o- c R>o"> a t the sites R. of the non- 
magnetic defects. 

In the large— N approac h 118 * 19 !, is represented by 
pseudofermion operators /J CT , f ra whose occupation num- 
ber n/(r) = J2<r(fLfrcr) obey the constraint n/(r) = 
N/2, with AT = 2 for S = 1/2. By adding the term 
Er« £ /( r )/r ) J r '° t° the Hamiltonian, the constraint 
(n/(r)) — N/ 2 can be enforced through the on-site en- 
ergy e f (r)2!122l The quartic interactions in Eq. (fil can 
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then be decoupled by introducing the mean fields 

= \ £tf"LO ; X(r,r') = ■ 

a a 

describing the local hybridization and the antiferromag- 
netic correlations, respectively. The resulting quadratic 
Hamiltonian can be diagonalized in real space [assuming 
periodic boundary conditions for an (L x L) lattice, where 
below, we take L = 41], allowing a self-consistent calcu- 
lation of s(r), x( r i r ')j an d £ /( r )- Below, we take for con- 
creteness J — 2t and I = 0.002i (unless otherwise stated) 
and thus study systems well inside the Kondo screened 
regime where s(r) ^ for all sites and magnetic fluctua- 
tions and the resulting c orrec tions to the mean-field level 
are expected to be wealP*l2ll 

Results We begin by considering a system with 17 
vacancies^ (~ 1% concentration of defects) and present 
in Figs, [lja) and (b) the relative variation of the hy- 
bridization As(r) = (s(r) — sq)/so, where sq is the 
hybridization of a clean lattice, and the corresponding 
absolute value of its Fourier transform (FT), |As(q)|, 
respectively. As predicted earlier^ and recently con- 
firmed by STS experiments on UR^SiJ^, As(r) exhibits 



isotropic oscillations with wavelength A s 



X F /2 



[see Fig. [TJa)], arising from the Fermi surface of the un- 
hybridized conduction band [see Fig. Wc)]. To under- 
stand the general momentum dependence of |As(q)|, and 
the effects of random disorder, we present in Fig. fijd) 
|As(qO for a single vacancy. A comparison of |As(qJ| in 
Figs, fl'b) and (d) shows that a finite concentration of 
defects leads as expected to a less well-defined momen- 
tum structure of |As(q)|. Moreover, a detailed analysis 
shows that the momentum dependence of |As(q)|, and 
in particular its maxima at q — (tt,tt), q = (0, 7r) and 
q = (tt, 0), are predominantly determined by the four 
strong peaks in |As(r)| in the immediate vicinity of the 
vacancy [see inset of Fig. [IJd)] . In |As(q)|, these peaks 
completely overshadow the long-distance, Xp/2 oscilla- 
tions in As(r) (whose amplitude is much smaller) such 
that their expected signature in |As(q)| at q s — 2tt/X s ss 
1.26 is not clearly observed [see Fig. fl^b)]. Finally, we 
present in Fig. [lje) the relative variation of the magnetic 
bond variable A%(r, r') = (%(r, r') — Xo)/xo [plotted at 
r m = (r + r')/2] where xo is the magnetic susceptibility 
for a clean lattice. As before^, we find that A%(r, r') ex- 
hibits strong oscillations along the lattice diagonal, aris- 
ing from the large degree of nesting of the hybridized 
F9^. As a result, the Fourier transform Ax(q m ) of 
Ax(r m ) shown in Fig. [TVf ) exhibits the anisotropic form 
of the hybridized FS [see Fig. Wlc)] (note that Ax(r m ) 
possesses lattice constant of <Xo/2, implying that for its 
first Brillouin zone, one has — 27r/ao < Qm' y ^ < 27r/ao). 

We next consider the differential conductance dl/dV 
and the absolute value of its FT, the QPI intensity | AT(q) | 
(for a derivation, see RefP^). It was previously showifl3l 
that the energy-dependent dl/dV lineshape sensitively 
depends on the ratio t / /t c (which is of the order of a few 
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FIG. 1: (color online) Kondo lattice with 1% of vacancies. 
Contour plot of (a) As(r), and (b) the absolute value of its 
Fourier Transform, |As(q)|. (c) Large and small Fermi sur- 
faces for an unperturbed lattice, (d) |As(q)| and As(r) (see 
inset) for a single vacancy. Contour plot of (e) A^(r m ) and 
(f) |Ax(q)|. 



percent.), where t c and t c are the amplitudes for tunnel- 
ing of electrons from the STS tip into the conduction and 
/-electron bands, respectively. While the overall mag- 
nitude of |A(q)| varies with tf/t c , we find that its mo- 
mentum dependence is rather insensitive to tf/t c , and wc 
therefore take for concreteness tf/t c = 0.03. In Figs. ^ a) 
and (b) we present the QPI spectrum, obtained for the 
system shown in Fig. [I] at two different energies. The 
presence of randomly distributed defects smears out the 
QPI spectra, as follows from a comparison with those 
spectra obtained for a single defect, shown in Figs. J2^c) 
and (d). At the Fermi level, E = 0, the QPI spectrum 
exhibits a peak at qo = (tt, tt) [see Fig.J^c)], arising from 
scattering between the almost parallel portions of the FS, 
as shown in Fig. [2je) . With decreasing energy, new peaks 
appear in the QPI spectrum away from the diagonal [q% 
in Figs. [2jb) and (d)] whose spectral weight quickly be- 
comes larger than that of the peaks, denoted by qi, close 
to the diagonal [Fig. [2je)]. This shift is a direct con- 
sequence of the spatial variation of the hybridization, as 
follows from a comparison with a QPI spectrum obtained 
within the Born approximation^! where the hybridiza- 
tion is spatially constant, and the main spectral weight 



FIG. 2: (color online) Contour plot of the QPI intensity 
iV(q)| for (a),(b) the system of Fig. [I] and (c),(d) a sys- 
tem with a single defect, for E — and E — — O.OOlt. (e) 
Scattering vectors dominating the QPI intensity and equal 
energy contours, (f) QPI intensity obtained using the Born 
approximation. 



FIG. 3: (color online) Kondo lattice with 1% of non-magnetic 
impurities with U = — Q.7t. Contour plot of (a) As(r), and 
(b) the absolute value of its Fourier Transform, |As(q)|. QPI 
intensity for a system with 1% defects (c),(d) and with a single 
defect (e),(f). 



still resides with the peaks (at qi) close to the diagonal 
[see Fig. [2jf )] . We thus conclude that the spatial oscilla- 
tions in the hybridization affect the QPI spectra, and in 
particular their spectral weight distribution, which is of 
great importance for the interpretation of experimentally 
obtained spectra, and the extraction of the underlying 
electronic band structure. 

The above results change qualitatively when one con- 
siders a system with 1% of non-magnetic defects [see 
Fig. [3] with scattering strength U — —0.7t, where the 
defects are located at the same positions as the vacancies 
of Fig. [IJ A comparison of As(r) shown in Fig. [3|a) with 
that in Fig. [ija) demonstrates that the introduction of a 
non-magnetic scattering potential significantly alters the 
spatial pattern of hybridization oscillations. In particu- 
lar, the oscillations acquire a larger amplitude, become 
much more isotropic, and the maxima in As(r) in the im- 
mediate vicinity of the defect are rotated by ir/4. These 
changes are particularly apparent in the Fourier trans- 
form, As(q) shown in Fig.^b), which exhibits an almost 
isotropic pattern. These changes arises from a strong 
feedback effect of the charge density on the hybridiza- 
tion: the s-wave form of the non-magnetic scattering po- 
tential leads to almost isotropic spatial oscillations of the 



conduction electron charge density^ (not shown) , which 
are reflected in those of As(r). As a result, the oscilla- 
tions in As(q) are dominated by 2fcj?-oscillations of the 
unhybridized Fermi surface [see black line in Fig. [3^b)] 
reflecting the strongly coupled nature of the system. The 
fact that the maximum intensity in As(q) is located at 
slightly larger momenta than 2kp arises from the expo- 
nential envelope of the spatial oscillations in Aa^ anc i 
the short decay length (£ « 2.2). Thus, while both 
vacancies and non-magnetic defects lead to spatial os- 
cillations in As with wavelength \ s = Xp/2, these os- 
cillations only become visible in the Fourier transform 
As(q) if their amplitude is enhanced by oscillations in 
the conduction electron charge density induced by the 
non-magnetic defects. These conclusions provide insight 
into the microscopic nature of defects: in particular, the 
recent observation of hybridization oscillations in 1% Th- 
doped URu2SiJ^, where the same 2fcj?-oscillations can be 
found in As(r) and As(q) suggests that the Th-atoms 
exert an appreciable non-magnetic scattering potential. 
The changes in As(q) are also reflected in the QPI spec- 
tra shown in Figs. p|c) and (d) [for comparison, we also 
present in Figs. [3^e; and (f) the QPI spectra for a system 
with a single non-magnetic defect] , which now differ sig- 
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FIG. 4: (color online) Kondo lattice with 1% of vacancies. 
Contour plot of (a) As(r), and (b) |As(q)| at T = 0.00002*. 
(c) As(r), and (d) |As(q)| at T = 0.00022*. (e) Equal energy 
contours at ifesTi^- (f) QPI spectrum for E — — 0.001* at 
T = 0.00022*. 



nificantly from those obtained in a system with vacancies 
[see Fig. [21. In particular, new peaks emerge in the QPI 
spectrum for E = [see black line in Fig.j3^e)] which re- 
flect those peaks found in iV(q) [see Fig. |3jb)] . We note 
that in the area between two (closely) placed defects [see 
yellow arrows in Fig. [3ja)] the spatial hybridization pat- 
tern is strongly affected by non-linear quantum interfer- 
ence effects. In particular, comparing the local hybridiza- 
tion pattern with that obtained by simply superposing 
the hybridization patterns of single, non-interfering de- 
fects, we find that these non-linearities in general sup- 
press large oscillations in the hybridization, i.e., in As(r). 
Finally, we verified that the above results remain qual- 
itatively unchanged for defect concentrations up to 2% 
and different spatial disorder realizations. 

Next, we demonstrate that even well below the Kondo 
temperature of the system, the spatial hybridization pat- 
tern and the QPI spectra can exhibit a significant tem- 
perature dependence. To this end, we again consider 
the system of Fig. [I] with 1% of vacancies but with a 
smaller J = 1.47* and I = 0.0002*. In Figs, ga) and 
(b) we present the resulting As(r) and As(q), respec- 
tively, at T\ = 0.00002*, which are similar to the results 
shown in Fig. [T] However, upon raising the tempera- 
ture to Ti — 0.00022*, we find a qualitative change in 



the spatial form of As(r) [see Fig. gc)]: the spatial ex- 
tent of the oscillations as well as their overall amplitude 
increases significantly. As a result, there also occurs a 
significant redistribution of spectral weight in As(q) [see 
Fig.gd)], which is now dominated by four peaks along 
the bond directions. The enhanced oscillations in the hy- 
bridization are now directly reflected in the QPI spectra, 
as shown in Fig. gf ) where we present the QPI pattern 
for E = — 0.001* which exhibits a very similar struc- 
ture to As(q). We note that these drastic changes in 
the hybridization and the QPI spectra occur well below 
the Kondo temperature Tr-, which we estimated from 
the vanishing of the hybridization for a clean system as 
Tjc « 0.004*, such that T x <T 2 «.T K . To understand 
these significant changes, it is necessary to consider the 
states in the Brillouin zone (BZ), which are excited at 
Ti.2, and which (roughly) lie between the equal energy 
contours E = ife^Ti^ shown in Fig. ge). At T\, states 
in only a small region of the BZ are excited [between the 
blue lines in Fig. pjfe)], leading to the pattern of As(r) 
shown in Fig. gaj! However, at T%, states in a much 
larger portion of the BZ are excited [between the red 
lines], leading to the significant changes in As and the 
QPI spectra shown in Figs.[4jc), (d) and (f), respectively. 
The strong temperature dependence discussed here is 
mainly an effect of the very weak /-electron dispersion, 
resulting from a small value of /, and will therefore de- 
crease with increasing strength of the antiferromagnetic 
interactions. 

Our results discussed above possess two important ex- 
perimental implications. First, it is apparent from the 
contour plots of As in Figs. [3^a) and |4jc) that already 
a small concentration of defects, indeed as small as 1%, 
can essentially disorder the hybridization in the entire 
system. To quantify this, we consider a site, r, as dis- 
ordered when the hybridization s(r) deviates by more 
than 1% from its value, so, in the clean system. We then 
find that for the case of non-magnetic defects shown in 
Fig. ga), 57.9% of the sites are disordered, whereas for 
the case shown in Fig. gc), 78.2% of the sites are disor- 
dered. This result provides an explanation for the strong 
disorder effects observed by Hamidian e* aZP^in (weakly) 
1% Th-doped URu2Si2. Second, while the hybridiza- 
tion can be strongly disordered even for small concen- 
tration of defects, as discussed above, the spatially aver- 
aged density of states at lu — 0, and thus also the specific 
heat of the system, is hardly affected. For example, for 
the two most disordered cases shown in Figs. [3ja) and 
gc), the spatially averaged density of states decreases 
only by about 8% for the system shown in Fig. 3]ja), 
while the change is less than 0.1% for the case shown 
in Fig. |4^c). These changes are quantitatively consistent 
with the changes seen experimentally in the specific heat 
of heavy fermion materials with defect concentration of 
19cP! Thus, we conclude that while the electronic struc- 
ture of heavy fermion materials can be heavily disordered 
already by small concentration of defects, the disorder's 
main signature appears in the hybridization and dl / dV 
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patterns, while the spatially averaged specific heat under- 
goes only modest changes, thus explaining the apparent 
contradiction between spectroscopy and thermodynam- 
ics. 
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